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ABSTRACT 


Results  are  presented  to  describe  the  laminar  flow  patterns  around  a  prolate 
spheroid  at  angles  of  attack  of  1,  2,  3,  and  30  degrees  and  complement  those 
obtained  previously  at  6  degrees.  They  were  obtained  by  solving  three- 
dimensional  boundary-layer  equations  with  a  combination  of  standard  and 
characteristic  box  methods  and  with  a  stability  criterion  to  ensure  numerical 
accuracy.  Emphasis  Is  placed  on  the  nature  of  separation  which,  1n  agreement 
with  experiment  but  contrary  to  some  theoretical  claims,  Is  shown  to  be  open 
for  all  angles  of  attack  and  to  be  coincident  with  a  particular  skin  friction 
line. 
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1.0  INTRODUCTION 


Prediction  of  separation  patterns  on  three-dimensional  bodies  is  Important  in 
engineering  practice  and  a  complete  solution  to  the  problem  requires  either  the 
solution  of  the  Navier-Stokes  equations  or  of  inviscid  flow  and  boundary-layer 
equations  with  Interaction  between  the  two.  The  calculations  are  difficult  and 
remain  to  be  performed  satisfactorily.  In  the  meantime,  attention  has  been 
directed  to  the  solution  of  model  problems  where  the  external  velocity  distri¬ 
bution  has  been  determined  from  an  analytic  expression  based  on  a  simple  body 
shape  with  three-dimensional  boundary-layer  equations  solved  numerically.  It 
should  not  be  expected  that  the  measurements  of,  for  example,  Han  and  Patel 
(1975),  Meier  and  Kreplin  (1980),  Meier,  Kreplin  and  Vollmers  (1983),  and 
Costis,  Polen,  Hoang  and  Telionis  (1988),  will  agree  in  detail  with  the  solu¬ 
tions  of  the  model  problem  but  the  latter  allows  the  development  of  essential 
numerical  techniques  which  pave  the  way  for  the  more  complete  solution  referred 
to  above.  It  can  be  expected  that  the  model  problem  will  emphasize  the  sensi¬ 
tivity  of  calculation  methods  in  the  region  of  separation  and  allow  the  nature 
of  flow  separation  to  be  identified;  it  is  this  aspect  which  is  the  major  con¬ 
sideration  of  this  report.  Reference  should  be  made  to  the  contribution  of 
Williams  (1977)  for  a  discussion  of  the  various  definitions  of  boundary-layer 
separation  in  three-dimensional  flows. 

An  important  aspect  of  the  subject  was  considered  by  Lighthill  (1963)  who  dis¬ 
cussed  the  nature  of  separation  in  three-dimensional  steady  laminar  flows.  He 
postulated  that  three-dimensional  separation  corresponds  to  a  skin-friction 
line  which  is  at  odds  with  the  earlier  suggestion  of  Maskell  (1955)  who  pro¬ 
posed  that  separation  is  an  envelope  of  limiting  streamlines.  The  argument  has 
continued  in  a  series  of  reports  where,  for  example,  Wang  (1975)  and  Cebeci  et 
al.  (1981)  have  sided  with  Maskell  and  more  recently  Cebeci  and  Su  (1988) 
reported  calculations  which  are  consistent  with  Lighthill 's  view.  These  recent 
calculations  were  performed  with  a  numerical  scheme  which  involved  systematic 
adjustment  of  the  grid  according  to  a  numerical  stability  criterion  which  had 
been  subjected  to  careful  evaluation  in  relation  to  an  extensive  range  of 
steady  and  unsteady  flow  problems.  They  corresponded  to  the  laminar  flow 
around  a  prolate  spheroid  of  thickness  t  »  1/4  at  6  degrees  angle  of  attack 
and,  as  shown  in  Fig.  1,  there  is  a  region  of  positive  crossflow,  followed  by 
a  substantial  region  of  negative  crossflow,  a  separation  line  and  two  terminal 
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Figure  1.  Flow  regions  for  a  =  6*.  according  to  Cebeci  and  Su  (1988). 


lines  beyond  which  solutions  of  the  boundary-layer  equations  could  not  be 
obtained.  Of  particular  note  was  the  separation  line  which  did  not  extend 
from  the  leeward  to  the  windward  side  Implying  open  separation  and  Is  In 
contrast  to  the  earlier  findings  of  Wang  (1975)  and  Cebed  et  al.  (1981). 

The  nature  of  the  problem  may  be  explained  In  terms  of  the  comprehensive 
examination  of  the  influence  of  angle  of  attack  by  Wang  (1976)  who  presented  a 
figure  which  showed  the  sequence  of  separation  patterns  on  the  prolate  spheroid 
for  a  range  of  angles  of  attack.  It  showed  that  (F1g.  2a)  the  flow  separates 
at  a  fixed  parallel  for  zero  Incidence,  that  a  small  Increase  of  Incidence 
(o  a  3")  tilts  the  separation  line  slightly  and  that  separation  Is  preceded 
by  a  small  region  of  weak  reversal  of  the  circumferential  flow  (Fig.  2b).  As 
the  Incidence  angle  further  Increases  (a  »  6*.  Fig.  2c),  the  leeside  separa¬ 
tion  point  moves  rearward  rather  than  forward  so  that  the  separation  line  Is 
bent  as  shown  1n  an  enlarged  view  (Fig.  3a).  The  upper  branch  AC  Inclines 
rearward,  whereas  the  lower  branch  CO  Inclines  forward.  Meanwhile,  the  rever¬ 
sal  of  the  circumferential  flow  starts  ahead  of  the  midbody  on  the  leeside 
(see  Fig.  2c).  At  higher  Incidence,  the  separation  line  breaks  so  that  the 
lower  branch  CO  extends  forward  to  become  an  open  separation  line,  as  shown  In 
Fig.  2d  and  enlarged  1n  Fig.  3b.  The  change  from  a  closed  separation  at  low 
Incidence  to  an  open  separation  at  higher  Incidence  thus  evolves  gradually  and 
It  Is  difficult  to  Identify  a  particular  Incidence  at  which  this  change  takes 
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—  SEPARATION  LINE 

--  CIRCUMFERENTIAL  -  FLOW  REVERSAL  LINE 
LIMITING  STREAMLINE 


(b)  o  -  3* 


Figure  2.  Sequence  of  flow  separation  on  a  prolate  spheroid,  according  to  Wang 
(1976). 


Figure  3.  Initiation  of  open  separation  on  a  prolate  spheroid  a  a  »  6®, 
according  to  Wang  (1976). 

place.  Even  for  6*  Incidence  (Fig.  2c  and  Fig.  3a),  It  Is  difficult  to  say 
whether  the  separation  line  Is  closed  In  the  section  BC. 

The  description  of  the  previous  paragraph  Is  In  some  contrast  to  the  measure¬ 
ments  of  Heler  et  a1.  (1983)  and,  even  though  exact  concurrence  with  calcula¬ 
tions  based  on  specified  pressure  distributions  Is  not  possible,  agreement  with 
general  features  of  the  flow  might  be  expected.  Figure  4  shows  the  separation 
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figure  4,  Separation  lines  derived  from  oil  flow  patterns  of  laminar  boundary 
layers,  according  to  Meier  et  al,  (1983). 

line  detected  by  Meier  et  al.  from  their  observation  of  the  bifurcation  of  oil 
flow  patterns  and  It  1s  evident  that  the  separation  Is  open  for  all  angles  of 
attack  greater  than  zero  and  in  general  agreement  with  the  calculations  of 
Cebed  and  Su  at  6  degrees.  A  similar  conclusion  was  reached  by  Costls  et  al. 
who  used  dyes  and  particles  to  visualize  the  flow  close  to  the  surface  of  a 
prolate  spheroid  of  fineness  ratio  four,  rather  than  a  ratio  of  six  as  In  the 
results  of  Wang,  Meier  et  al.  and  Cebed  and  Su.  The  visualization  patterns 
show  that  their  separation  was  open  for  angles  of  attack  from  6  to  30  degrees 
and  less  readily  defined  at  3  degrees. 

Thus  the  purpose  of  the  present  report  Is  to  present  calculations  performed 
with  the  procedure  used  by  Cebed  and  Su  (1988)  and  to  determine  If  the  dif¬ 
ferences  between  their  results  at  a  =  6*  and  the  description  of  the  previous 
paragraph  also  exist  at  other  angles  of  Incidence.  At  the  same  time,  the 
results  are  examined  In  terms  of  the  postulations  of  Lighthlll  and  Maskell  to 
determine  the  nature  of  the  separation  line.  To  achieve  these  purposes, 
results  are  presented  for  Incidence  angles  of  1,  2,  3,  and  30  degrees. 
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The  following  section  provides  a  brief  description  of  the  method  used  to  solve 
the  three-dimensional  boundary-layer  equations  and  emphasizes  the  special  pre¬ 
cautions  taken  In  the  vicinity  of  separation.  The  results  are  presented  In 
Section  3,  which  discusses  them  1n  the  context  of  the  previous  studies  referred 
to  above.  The  report  ends  with  a  summary  of  the  more  Important  conclusions. 
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2.0  CALCULATION  METHOD 


The  calculation  method  used  by  Cebed  and  Su  (1988)  considers  the  conservation 
equations  for  mass  and  for  x-  and  9-momentum  In  curvilinear  orthogonal 
coordinates  and  obtains  their  solution  by  a  finite-difference  method  based  on 
Keller's  box  scheme  for  a  prolate  spheroid  defined  by 


subject  to  an  external  velocity  distribution  given  by  the  Inviscid  flow  theory. 
To  avoid  the  singularity  at  the  nose  E  (=  x/a)  =  -1 ,  In  a  region  -1  £  E  £  E^, 
they  used  the  procedure  of  Cebed  et  al.  (1980)  which  employs  transformations 
that  allow  the  boundary-layer  equations  to  be  solved  without  numerical  diffi¬ 
culties.  For  E  >  Eq,  they  used  two  different  versions  of  the  box  scheme 
depending  on  the  complexity  of  the  flowfleld.  In  regions  of  positive  crossflow 
velocity,  they  used  the  standard-box  scheme  and  the  characteristic  box  scheme 
In  regions  of  negative  crossflow  velocity. 

The  accuracy  of  the  calculations  performed  with  the  standard-box  method  Is  well 
established  and  that  of  the  characteristic-box  method  was  examined  by  Cebed 
and  Su  for  regions  of  positive  crossflow  and  was  shown  to  be  equally  satis¬ 
factory.  In  regions  of  negative  crossflow  and  In  the  vicinity  of  the  separa¬ 
tion  line,  however,  tne  results  depended  upon  the  finite-difference  grid  In  a 
way  which  necessitated  an  Intelligent  procedure  for  grid  spacing.  To  discuss 
this  point  further,  let  us  consider  the  net  shown  In  Fig.  5  at  a  distance  y 
from  the  surface  and  assume  that  the  solutions  with  the  characteristic  box 
scheme  originate  on  the  leeward  line  of  symmetry  and  march  towards  the  wind¬ 
ward  line  of  symmetry.  The  symbol  x  denotes  the  location  where  the  solution 
Is  Icnown  and  the  symbol  o  denotes  the  location  where  the  solution  Is  to  be 
found.  The  backward  characteristic  from  point  P  Is  1n  the  local  streamline 
direction  and  Intersects  the  E^  ^  Une  at  E  when  there  Is  a  positive  cross- 
flow  velocity  and  at  F  when  the  crossflow  velocity  Is  negative.  Since  the 
characteristic  box  computes  the  region  EFP,  which  Is  known  as  the  domain  of 
dependence  of  point  P,  It  ensures  that  the  domain  of  stable  computations  can 
be  determined  a  priori  and  this  Is  achieved  by  determining  the  ratio  P(rA0^/A0) 
which  Is  also  equal  to  (w/uKh^/h^XdE/AQ)  and  requiring  that  It  remains  small 
during  the  calculations:  Cebed  and  Su  (1988)  refer  to  this  requirement  as  the 
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X  KNOWN 
O  UNKNOWN 


Figure  5.  Finite-difference  notation  for  the  characteristic-box  scheme:  o, 
unknowns,  x  known. 


stability  criterion  of  the  characteristic  box  scheme  since  It  characterizes 
the  flow  angle  and,  therefore,  the  grid-interval  required  to  ensure  that  the 
zones  of  dependence  and  Influences  are  properly  represented. 
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3.0  RESULTS 


As  1n  the  study  of  Cebeci  and  Su  (1988),  the  calculations  for  a  given  angle  of 
attack  were  performed  by  identifying  three  separate  regions  on  the  body.  In 
the  region  upstream  of  the  line  of  zero  c^  (Region  A),  all  u-  and  w-velocities 
are  positive;  in  the  region  between  this  l?ne  and  the  two  terminal  lines  1  and 
2  (Region  B).  u  is  positive  and  w  is  negative  near  the  surface;  and  in  the 
region  downstream  and  between  the  two  terminal  lines  (Region  C)  u  and  w  are 
negative  near  the  surface  and  positive  away  from  the  surface.  The  numerical 
calculations  in  Region  A  made  use  of  the  regular  box  scheme  with  comparatively 
large  step  sizes  in  the  streamwise  direction  E  and  those  in  Region  B  made 
use  of  the  characteristic  box  scheme  with  a  great  deal  of  attention  was  paid 
to  the  choice  of  the  net  in  the  circumferential  and  streamwise  directions. 
Again,  three  subregions  were  identified  to  perform  the  calculations.  The 
first  subregion  started  at  where  the  w-velocity  became  negative  and 
extended  to  E^  where  the  streamwise  wall  shear  (du/dy)  first  vanished.  The 

r  W 

second  and  third  subregions  began  at  the  first  separation  point  Ep  and  were 
bounded  by  the  leeward  line  of  symmetry,  the  terminal  line  1,  the  windward  line 
of  symmetry  and  terminal  line  2. 

Following  the  routine  calculations  in  Region  A,  the  calculations  in  Region  B 
were  started  on  the  windward  line  of  symmetry  and  continued  towards  the  leeward 
line  of  symmetry  with  the  standard  box  scheme  where  w  is  positive  and  the 
characteristic  box  scheme  where  w  is  negative.  The  step  lengths  in  the  circum¬ 
ferential  and  streamwise  directions  were  chosen  to  satisfy  the  stability  param¬ 
eter  0.  A  typical  streamwise  grid  for  uniform  circumferential  grid  of  AS  »  2.5 
is  shown  in  Table  1  for  a  »  1,  2  and  3  degrees. 

The  previous  study  of  Cebeci  and  Su  (1988)  for  a  *  6*  showed  that  the  choice 
of  step  lengths  in  region  6  was  very  important  to  the  accuracy  and  the  stabil¬ 
ity  of  the  solutions.  If  they  were  not  chosen  properly,  the  solutions  oscil¬ 
lated  and  sometimes  led  to  breakdown  of  the  calculations.  The  severity  of  the 
oscillations  depended  on  the  magnitude  of  the  flow  reversal  and,  when  it  was 
small,  so  were  the  oscillations.  In  regions  close  to  the  flow  separation  line, 
however,  the  solutions  became  more  sensitive  to  the  criterion  set  by  the 
B-parameter  and  required  that  the  value  of  0  remained  relatively  small.  A 
precise  identification  of  the  first  (5p,0p)  location  on  the  body  where 
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Table  1 

Axial  step  lengths  for  &6 

=  2.5  degrees 

a  s 

r 

a  =  2 

• 

o  = 

3* 

i 

AE  X  10 

^  I 

AE  X  10 

1  i 

AE  X  1i 

0.0000-0.51 

1.0 

0.0000-0.430 

1  .0 

O.aOO-0.430 

1 .0 

0.5100-0.60 

0.5 

0.4300-0.550 

0.5 

0.430-0.500 

0.5 

0.6000-0.61 

0.25 

0.5500-0.590 

0.25 

0.500-0.545 

0.25 

0.6100-0.62 

0.125 

0.5900-0.606 

0.10 

0.545-0.560 

0.10 

0.6200-0.625 

0.1 

0.6060-0.606 

0.05 

0.560-0.572 

0.05 

0.6250-0.630 

0.05 

0.6080-0.6095 

0.025 

0.572-0.586 

0.025 

0.6300-0.6355 

0.025 

0.6095-0.6110 

0.01 

0.586-0.611 

0.05 

0.6355-0.6380 

0.01 

0.6110-0.6120 

0.025 

0.611-0.640 

0.1 

0.6380-0.6410 

0.025 

0.6120-0.6140 

0.05 

E  >  0.640 

0.125 

0.6410-0.6440 

0.05 

0.6140-0.630 

0.10 

0.6440-0.6500 

0.1 

E  >  0.630 

0.25 

0.6500-0.6550 

0.125 

1  >  0.6550 

0.25 

the  streamwise 

« 

wall  shear  f  vanished  and 

N 

where  g 

was  negative  was 

very  Important 

In  the 

calculations  since 

that  point  corresponded  to 

the 

definUlon  of  the  open  flow  separation  line  on  the  prolate  spheroid. 

In  view  of  the  experience  gained  In  computing  o  .  6*  case,  our  studies  were 
conducted  with  step  lengths  chosen  to  satisfy  the  stability  parameter  B  (see 
Table  1)  and  the  solutions  were  examined  to  ensure  that  they  were  smooth  and 
free  of  oscillations.  The  grids  used  to  compute  the  flow  at  a  >  1 ,  2  and  3 
degrees  were  essentially  the  same.  For  this  reason,  the  details  of  the 
procedure  are  discussed  for  a  »  3  and  the  remaining  two  will  be  summarized. 

For  o  =  3*.  the  flow  reversal  1n  the  circumferential  flow  direction  takes 
place  at  1.  =  0.08  and  the  boundary  of  the  line  c,  =0  varies  with  I  as  shown 

A  I  A 

1n  Fig.  6.  The  first  ^-location  where  the  streamwTse  wall  shear  parameter 
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Figure  6.  Flow  regions  for  o  «  3*. 

II 

f  vanishes  occurs  at  |  =  0.57375,  0p  »  100*.  Figure  7  shows  the  variation 

of  f^  with  I  and  0,  for  several  values  of  0  and  and  we  note  that  f^ 

vanishes  at  5  =  0.765  on  the  windward  line  of  symmetry  (0  «  0),  and  moves 

N 

upstream  with  Increasing  0  until  100*  where  I  *  0.57375.  For  values  of  0  greater 

M 

than  0  a  100*.  the  ^-location  where  f^  Is  zero  Increases  to  5  »  0.595  for  0 
=  140*  and  to  5  »  0.617  for  0  »  180*.  The  variation  of  f^  with  0  for  given 
^-locations  on  the  body  Is  shown  1n  Fig.  7b  which  Indicates  that,  for  values  of 
5  <  0.57375,  the  strearowlse  wall  shear  decreases  monotonically  with  Increas¬ 
ing  0  up  to  a  certain  ^-location  after  which  It  begins  to  Increase.  A 
sharp  decrease  In  Its  magnitude  takes  place  at  5  »  0.57375  for  0  »  100*  where, 
after  a  steep  decrease  1n  f  ,  It  shows  a  rather  sharp  Increase.  We  note  from 

H  ^ 

the  behavior  of  f  at  this  angle  of  attack,  and  those  at  lower  angles,  that 
this  rapid  decrease  of  f^  at  a  certain  %  and  0  values  Is  much  less  than  those 
computed  at  higher  angles  of  attack,  I.e.  o  •  6*.  As  a  result,  the  calcula¬ 
tions  were  easier  to  perform  at  the  lower  angles. 

We  see  from  Fig.  7,  that  the  computed  wall  shear  values  are  smooth  and  finite 
and  free  of  oscillations.  The  results  of  Fig.  8  support  this  conclusion  In 
that  the  magnitude  of  the  stability  parameter  8  Is  shown  to  be  small  In  the 
circumferential  direction  at  various  ^-locations  on  the  body. 
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hear  fy 


with  (a)  (  and  (b) 


respectively,  for  o  »  3*. 


Figure  8.  Variation  of  the  stability  parameter  3  with  0  at  general  values 
of  5  for  a  =  3*. 

Figure  9  shows  the  approximate  direction  field  of  the  limiting  streamlines 

H  H 

Qw^^w  indicated  by  dashed  lines)  and  the  Integrated  trajectories  (Indicated 
by  solid  lines)  In  the  0-5  plane.  Figure  10  shows  the  limiting  stream- 

H  H 

lines  obtained  by  Integrating  the  direction  field  of  the  wall  shear 
as  a  function  of  5  and  0  plotted  on  the  body.  We  note  that  the  limiting  stream¬ 
lines  which  were  parallel  to  the  ^-direction  at  o  «  0  are  now  deflected  from 
the  5-dlrectlon  due  to  the  effect  of  the  pressure  gradient  In  the  0-dlrectlon. 

As  In  the  case  of  a  »  6*,  the  calculations  lead  to  definitions  of  terminal 
lines  1  and  2  which  correspond  to  the  breakdown  of  the  calculations  due  to  the 
singular  nature  of  the  boundary-layer  equations.  In  these  regions,  the  slopes 
of  the  limiting  streamlines  change  rapidly  and  approach  *90*  as  the  solutions 
break  down.  It  Is  clear  from  this  figure  and  from  Fig.  6,  that  the  first 

H 

point  at  5  »  0.57375,  0  ■  100",  which  corresponds  to  negative  g^  and  zero 

M  w 

f^  can  be  regarded  (as  In  the  case  of  a  »  6*)  as  the  beginning  of  the  separa¬ 
tion  line  for  o  ■  3*  and  can  be  used  to  Identify  the  skin-f rictlon  line  as  the 
separation  line  In  three-dimensional  flows.  The  behavior  of  the  sk1n-fr1ct1on 
lines  on  both  sides  of  the  separation  line  Is  again  similar  to  that  for  a  >  6* 
discussed  by  Cebeci  and  Su  (1988). 

Figures  11  to  13  show  the  results  for  a  *  2*  and  1*  and  the  first  flow 
reversal  In  the  circumferential  flow  direction  takes  place  at  *  0.17  and 
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Figure  9.  Approximate  direction  field  of  the  limiting  streamlines  (dashed 

lines)  and  Integrated  trajectories  (solid  lines).  The  long  dashed 
line  SP  denotes  terminal  line  1,  PN  terminal  line  2,  and  the  solid 
line  PN  the  separation  line. 


Figure  10.  Behavior  of  the  limiting  streamlines  on  the  body  fo 
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Figure  11.  Flow  regions  for  a  »  2*.  5$  ^nd  ^|i|  correspond  to  the  separation 
locations  on  the  leeward  and  windward  sides,  respectively. 


TERMINAL 
LINE.  T^ 


SEPARATION 

UNE 


Figure  12.  Flow  regions  for  a  »  I*.  5$  af**!  5n  correspond  to  the  separation 
locations  on  the  leeward  and  windward  sides,  respectively. 


1 


Figure  13.  Variation  of  the  streamwise  wall  shear  fy  with  I  at  different 
values  of  9  for  a  *  1*. 


at  =  0.290,  respectively.  The  first  ^-location  where  f^  vanishes  occurs  at 
=  0.6098,  Op  =  130“  for  a  .  2*,  and  at  =  0.6374,  Op  .  175*  for  a  .  1*.  As 
for  a  =  3*,  the  figures  define  two  terminal  lines  and  a  separation  line  which 
Is  open  and  show  that  the  beginning  of  the  open  separation  line  approaches  the 
separation  location  on  the  leeward  line  of  symmetry  with  Increasing  angle  with 
open  separation  closed  at  a  »  0.  Figure  13  shows  the  variation  of  the  stream- 

M 

wise  wall  shear  f^  with  0  at  various  ^-locations  on  the  body  for  o  =  1*  and 
that,  as  for  a  »  3*  and  the  solutions  are  free  of  oscillations. 

Figure  14  shows  the  open  separation  and  c^  *  0  lines  where  the  first  flow 
reversal  In  the  circumferential  direction  ?akes  place  at  5  »  -0.92  on  the 

41 

leeward  plane  of  symmetry,  and  where  f^  vanishes  at  5  *  -0.8453  at  6  »  141.5* 
for  a  »  30*.  The  calculations  for  this  angle  of  attack  were  more  difficult 

than  those  for  o  »  6*  and  It  was  necessary  to  use  a  nonuniform  grid  In  both 
5  and  0-d1rect1ons,  as  shown  In  Tables  2  and  3.  It  Is  evident  from  Fig. 

15  that  the  magnitude  of  8  Is  relatively  small  for  values  of  5  <  -0.875 
and  begins  to  Increase  sharply  thereafter.  For  I  *  -0.8453,  the  magnitude 
of  8  reaches  0.57  at  0  *  141.5*  and  It  was  possible  calculate  up  to  5  » 

-0.50  In  the  second  subregion.  Overshoots  In  the  velocity  profiles  (both  In  u 
and  w)  began  to  occur  at  5  *  -0.50  and  became  progressively  worse.  Refining 
the  net  did  not  help,  and  the  calculations  became  very  difficult  to  perform 
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s 


-TERMINAL  LINE.  T, 


-TERMINAL  LINE.  Tj 


-SEPARATION  LINE  (VERY  CLOSE  TO  Tj) 


Figure  14.  Flow  regions  for  a 


AO  ISO  IM  170 

•  (d« 

Figure  15.  Variation  of  the  stability  parameter  0  with  0  at  several  values 
of  5  for  a  ■  30*. 


with  Increasing  I  where  the  circumferential  edge  velocity  became  bigger  and 

N 

the  boundary- layer  thicker.  Figure  16  shows  the  variations  of  f^  with  5  and 
9,  for  several  values  of  0  and  I,  respectively,  and  Indicates  trends  similar 

N 

to  those  observed  with  a  *  6*.  We  note  from  Fig.  16a  that  f^  vanishes  at 
3  0.980  on  the  windward  line  of  symmetry  (0  «  0),  and  at  f  >  0.948  on  the 


Table  2.  Axial  Step  Length  for  a  =  30 


i 

<- 

■0.9 

0.25 

-0.9000 

-  -0.8640 

0.10 

-0.8640 

-  -0.8520 

0.05 

-0.8520 

-  -0.8500 

0.04 

-0.8500 

-  -0.8470 

0.03 

-0.8470 

-  -0.8456 

0.02 

-0.8456 

-  -0.8446 

0.01 

-0.8446 

-  -0.8418 

0.02 

-0.8410 

-  -0.7000 

0.10 

>-0.7000 

0.20 

Table  3.  Circumferential  Step  Lengths  for  o  »  30* 


1  < 

-0.8418 

1  > 

-0.8418 

8 

8 

4® 

0* 

- 

30* 

5* 

0* 

- 

30* 

5* 

30* 

- 

100* 

2.5* 

30* 

- 

100* 

2.5 

O 

o 

o 

- 

110* 

2* 

100* 

- 

no* 

2* 

no* 

- 

139* 

1* 

no* 

- 

139* 

1* 

139* 

- 

144* 

0.5* 

139* 

- 

141* 

0.5 

144* 

- 

160* 

1* 

141* 

- 

170* 

3* 

160* 

- 

166* 

1.5* 

170* 

- 

180* 

5* 

166* 

- 

170* 

2* 

170* 

180* 

2.5* 
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leeward  side.  The  first  Interior  point  where  f  vanishes  occurs  at  5-  = 

-0.8453  at  0^  =  141.5“.  The  variation  of  f^  with  0  for  a  given  5-locatlon, 
see  Fig.  16b,  shows  that,  the  wall  shear  undergoes  a  drastic  change  near  the 
region  of  flow  separation,  decreasing  and  Increasing  sharply  at  0  =  141.5“ 
and  this  explains  the  need  to  employ  a  very  fine  grid  In  this  region.  In  both 

It 

figures,  the  behavior  of  f  Is  smooth  and  free  of  oscillations  and 

w 

Indicates  that  the  results  are  accurate  and  reliable. 

Figure  17  shows  the  approximate  direction  field  of  the  limiting  streamlines 

(Indicated  by  dashed  lines),  the  Integrated  streamlines  (Indicated  by  solid 

lines),  the  c^  =0  line,  separation  line  and  terminal  lines  1  and  2.  There  Is 

clearly  no  problem  In  computing  the  flow  In  the  third  subregion  after  flow 

separation  occurs  at  5  =  -0.8453  because  the  negative  cross-flow  region  Is 

very  small.  Consistent  with  the  results  for  a  »  1  to  3*  and  6“,  there  Is  a 

separation  line  between  the  two  lines  representing  c,  »  0  and  the  terminal 

*0 

line  1.  The  approximate  streamlines  were  obtained  by  Integrating  g^/f^  as  a 
function  of  5  and  0  with  Integration  starting  at  5  »  -0.860  and  their 
behavior  Is  similar  to  those  for  a  »  6*. 

The  separation  lines  of  Fig.  18  may  be  compared  with  those  of  Figs.  2  to  4  with 
the  clear  Impllcaton  that  the  present  results  agree  well  with  the  measurements 
of  Meier  et  al.  and.  In  the  main,  with  those  of  Costls  et  al.  The  separation 
Is  open  for  all  angles  of  attack  greater  than  zero,  at  least  for  the  prolate 
spheroid  of  fineness  ratio  six.  It  should  be  noted  that  the  side  view  pre¬ 
sented  on  Fig.  18  shows  a  small  distance  of  open  separation  for  angles  of 
attack  of  1 ,  2  and  3  degrees  so  that  great  care  Is  likely  to  be  required  In 
the  Interpretation  of  near-surface  visualization,  as  In  the  case  of  Costls  et 
al. 
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4.0  CONCLUOIN6  REMARKS 


The  results  of  the  preceeding  section  confirm  that  separation  Is  open  for  all 
angles  of  attack  considered  and  described  by  a  line  which  corresponds  to  a 
skin-friction  line,  as  suggested  by  Lighthlll.  The  lines  of  zero  c^  occur 
upstream  of  separation  and  terminal  lines  are  Identified  beyond  whlc^  the 
boundary- layer  equations  cannot  be  solved  with  a  prescribed  pressure  distri¬ 
bution.  It  Is  evident  that  great  care  Is  required  to  ensure  numerical  accuracy 
and  that  this  has  been  achieved  by  a  combination  of  the  standard  and  charac¬ 
teristic  box  methods  with  an  essential  stability  criterion. 
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